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Abstract 

We show how the partition function of a network of parallel superconducting wires weakly coupled 
together by the proximity effect, subjected a vector potential along the wires can be mapped onto 
N-distinguishable two dimensional quantum-mechanics problem with a perpendicular imaginary 
magnetic field. Then, we show, using a mean field approximation, that, for a given coupling, there 
is a critical temperature for onset of inter-wire phase coherence. The transition temperature Tc is 
plotted on both cases for non-magnetic and a magnetic field perpendicular to the wires. 

PACS numbers: 
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I. INTRODUCTION 



There has been considerable recent interest in thin wires that undergo transitions into an 
ordered state, such as superconducting or ferromagnetic. For example, a recent experiment 
[l| has suggested that single-walled carbon nanotubes (which have diameters of only about 4 
A) are superconducting up to temperatures as high as 20 K. Because these tubes are so thin, 
they behave very much like one-dimensional superconductors. It was therefore proposed |l| 
that they could be described by a complex order parameter ip{z) which varies only in one 
dimension, say the z direction, i.e. along the tube, iplz) might represent the complex energy 
gap, or, in a different normalization, it could represent the condensate wave function in a 
BCS superconductor. 

Moreover, there have been many experiments for investigating superconductivity on 
nanowires. Ropes of carbon nanotubes between superconductingelectrodes can show super- 
conductivity due to the proximity effect of the electrodes . Furthermore, supercon- 
ductivity on carbon nanowires connected to normal contacts, has been observed io). On 
the other hand, superconductivity of nanowires of Zn or Sn has been investigated [7|, [Sj. 

Fluctuations are, of course, especially important in one- dimensional systems. It was 

n 

shown many years ago by Scalapino et al [9] that classical fluctuations in one dimension 
could be treated exactly, within the context of a Ginzburg-Landau (GL) free energy func- 
tional. Their treatment involved mapping the GL functional onto a single-particle quantum 
mechanics problem, using an exact connection between the classical partition function and a 
path integral treatment of the quantum mechanics problem. These authors showed that clas- 
sical fluctuations could give rise to a non-zero order parameter even above the GL transition 
temperature. This mapping was extended to treat Josephson-coupled thin wires 
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111. 



However, in the mapping, the effect of a magnetic field was ignored. In the case of a non- 
zero perpendicular magnetic field, we show that the partition function for the wires maps 
onto a certain zero-temperature quantum mechanics problem in two dimensions with an 
effective imaginary perpendicular magnetic field, which brings to a non-Hermitian quantum 
mechanics problem. 



The non-Hermitian problem in physics has not been new recent 



cesses can be described by non-Hermitian Liouville operators 



12 



13 



y. Nonequilibrium pro- 



14| . The non-Hermitian 



quantum mechanics are well studied in order to study the pinning of magnetic flux lines in 
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high temperature superconductors 




The remainder of this paper is organized as follows. In Section [TTl we describe our formal- 
ism and mapping. In Section llllt we give our numerical results including phase diagrams. 
This is followed by a concluding discussion and an outline of possible future research. 



II. FORMALISM 



A. Mapping to a quantum mechanics problem for interacting superconducting 
wires when B is perpendicular to the wires 

Let us consider a network of N parallel superconducting wires in a non-zero vector poten- 
tial. We assume, for convenience, that these wires all have the same GL parameters, though 
the formalism can easily be generalized to the case when the parameters are different. Then 
the partition function can be written as a functional integral over the complex order 
parameters ipii^i), ■■■''Pn{.Zm): 

VMzi)-'DiJN{zN) exp{-(3F[ij^{zi), ...^n{zn)]}. (1) 

We assume that the free energy functional is the sum of two parts: a single-wire term Fg 
and a term describing inter- wire interactions, which we denote Fj^j. The single- wire term 
will just be the sum of general GL equation for each wire: 

N 

Fs = J2FGLmz^)]. (2) 



i=l 



Here, 



FGLmz^)] = [^1 (-V-^j^(;^)P + a|^(z)r + 7|^(^)r + ^S]ciz, (3) 

where a, 7, and m* are material-dependent (and possibly temperature-dependent) coeffi- 
cients. Commonly, it is assumed that 7 is positive and that a = a'{T — Tc), where T is 
the temperature, is the critical temperature, and a' is greater than zero. Also, S is the 
cross-sectional area of the sample, but for one-dimensional wire we may ignore this term. 
For the interaction term, we assume a form similar to that used by Lawrence and Doniach 

n 

for interacting superconducting layers |18l], namely 



int 



r^max 

K,m^)-Uz)\\ (4) 



where Zmax is the length of the wires. Basically, we are assuming that there is a Josephson 
coupling of strength Kij between different wires, but at the same point along the length, z. 
We choose a gauge such that the vector potential is parallel to the superconducting wires, has 
only z component and independent of z. When a wire is a loop, a vector potential is related to 
the total flux $ through the loop, = ^1 z^ax- In this case, using ipi{z) = ipiR^z) +iipii{z), 
Fs and Fint take the forms 



and 



K,,{\^i{z)\^ + \^j{z)\^ - 2(V'iijVii? + Aii^ji)), (6) 
where ip'{z) = dil){z) / dz. Finally, the partition function takes the form 

Z= lY,'^^iR'^'Pii''M-l3F[iJiR,M). (7) 

where we use ■j/'j = ipiR + iipu. 

We now show that eqs. ([5]), ([6]) and (JTj) for Z are actually equivalent to a quantum 
mechanical problem of a N distinguishable particles in N distinct quantum wells in two 
dimensions in the presence of a perpendicular magnetic field. In order to simplify our 
argument, we consider the case of single particle with mass, m and a charge e* subjected to 
a 2D potential, V{x,y). The density matrix of a two-dimensional system, using ip^ and ip^ 



are boundary condition at initial and final time, can be written 19 1 



(^^|e-^/^|^^) = C Vx{r)Vy{T)e^Y>{-\s[x{r),y{r)]}, (8) 
1,1 n 



where 



S = / + y") + Vix,y) - i-A ■ vjdr. (9) 

For the given B = Bej/z with the gage 

Aff = ^{xy-yx), (10) 

this S becomes 

rf^effH ™ p*B ff 

S= drl^ix" + y") + Vix,y) - ^-^{xy' - yx% (11) 



Q.M. 


S.C. 


T 


z 


Pi = {xi{u),yi{u)} 




E 


p 

^max 






Vi{xi,yi) 




m 




^eff 




J'ij 


IT 



TABLE I: Correspondence on the mapping between Q.M. and S.C. (in each case, the left-hand 
variable corresponds to the parameters on the quantum mechanics problem and the right hand 
variable corresponds to the parameters on the superconductor wires) 

This is a similar equation to the partition function of the superconducting wires. 

In order to simplify this mapping, we use the suitable dimensionless form. 
Co^'^i^iR, and -ipiy = ^o^'^ipu- Then we can make the identifications of Table [B 

We find that the magnetic field has two effects: (i) it determines an effective perpendic- 
ular magnetic field in which the equivalent quantum-mechanical particle moves; and (ii) it 
changes the quadratic part of the effective potential. The Hamiltonian for the analogous 
quantum problem is 

H = h^ (p- + ^fy) ' + i (p^y - ' + ^^(/^^)] + E 2^^>i - p. r> (12) 

i=l ^ ^ ^ ^ <ij> 

where and Piy are momentum operators of x and y components of ith particle, respec- 
tively. 

B. Probability distribution of the order parameter 

We consider the probability distribution of the order parameter, which corresponds to 
the probability distribution of particles in quantum mechanics. In order to simplify our 
discussion, we consider single wire case. The probability distribution function of the order 
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parameter can be defined as 

where Z = {ilj^\e^^^^\il)^) and IV'^) represents the boundary condition at r = and 
represents the boundary condition at r = L^-. Using the eigenstates of the Hamiltonian, 
H\n) = En\n), the probability can be written as 

PiPir)) = |$^(^^|m)(m|p(r))(p(r)|n)(n|^Oe-^(^^-^)e-T- (14) 

m,n 

with 

Z = J2{^^\n){n\ilj')e-^\ (15) 

n 

Explicitly, the expectation value of operator, p at the distance r from the bottom of the 
wires is given by 

{P)r = ^{^le-^^'-^^ J dp{rMr))p{pir)\e-^^m, (16) 

where p\p) = p\p). 

The case of periodic boundary condition, our problem can be simplified. If corresponds 
to and summed over all possible this configuration, the density matrix can be written as. 

m,n I 

= ^J2^Mpir)){p{r)\n)e-'^'^\ 

n 

where Z = ^j{'ip^\n){n\'i/j^)e~^^^ — e~"^^^. So, if the wire is actually in the form 
of a loop, which means the boundary conditions ■0(0) = ip{zmax), our problem corresponds to 
this statistical mechanics. Of course, in the limit of a very long wire, the periodic boundary 
condition imposed by the loop should become unimportant. 

In the case of the periodic boundary condition for single wire, we can see qualitative 
behavior of order parameter. The average gap in the GL problem (denoted A(t)) corresponds 
to the mean distance (p) in the quantum-mechanical problem, i.e. 

ip) - A(i), (17) 

where A(t) — A.(t)/^Q^'^. At much lower temperature than the critical temperature T^, the 
mean distance from the origin of the particle approaches the value predicted for the quantum 
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problem in the limit of infinite mass, i.e. the value of p for which the quartic potential is 
a minimum although the magnitude of these gaps at T = are different. The function 
a/1 — t, is the classical solution, i.e., in the case when thermal fluctuations in the GL case 
are negligible. These fluctuations do indeed become very small when T ^ 0, because in 
this regime, the effective potential rises steeply above its minimum, and the (p) becomes 
very close to the value that minimizes the GL free energy. When (p) has this value, the 
corresponding value for A(t) is 

A(t) = V^V + v^/ = ^^^v/r^ = Hmt). (18) 

where A(0) is the gap at T = 0. These considerations may suggest that we can approximate 
g{t) = A(t)/A(0) = Vl^t. 



C. Phase only model and mean-field approximation 

This system will undergo a phase transition into a phase-ordered state below a critical 
temperature Tc which is distinct from (and lower than) the single wire mean-fleld tran- 
sition temperature T^. To do this, we consider a simplifled, "phase-only" version of this 
Schrodinger equation (fT2!) . We assume that the magnitudes pi of the variables Xj are flxed 
at the values which minimize the single-wire GL free energy, i.e. pi = po ( |T8l) . All terms 
in the Hamiltonian involving d/dpi can be ignored in this phase-only model. The effective 
Hamiltonian f[T^ then becomes 

^ - - E ^ij - E - 2 5: - cos - (19) 

where the sum runs over distinct nearest neighbor pairs. This is the well-known quantum 
XY model, which exhibits a quantum phase transition at a critical value. 

The mean fleld approximation can be applied to this Hamiltonian, assuming that Jij = J 
for only nearest neighbors, by replacing the second term according to the prescription 

cos — (f)j) = 2 COS (f)i {cos (f)) — (cos</))^, (20) 
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where we are supposing (sin 0) = because of the symmetry. Thus, 

2 JPoi^ — cos {(pi — (f)j)) = 2 Jpl (l — 2 cos0i(cos0) + (cos 0)^) 

= — 4z„ Jpq(cos0) ^^cos0j + 2 ^Po(l + (cos0)^) 
= — 4z„Jpq(cos0) ^cos0j + 2znNJpl{l + (cos0)^), 

where Zn is the number of nearest neighbors in the lattice. Thus, the effective Hamihonian 
corresponding to eq. f|T9l) becomes a following Schrodinger equation: 



" " — 42;„p^ J (cos 0)cos 0i+ 22;„ Jp^ ( 1+ (cos 0) (0i ) = {(pi ) 

21) 



2mpQ 90^ 2mc i 



We consider the self consistent equation for cos0 on the periodic boundary condition. 
The mean field theory is defined by the self-consistency requirement on (cos0): 

/ n _ ^„e-/^-//^"(V>„(0,)|cos0#„(0,)) 

For example, when the wires are sufficiently long where only the ground state contribution 
may be important, the self-consistent condition becomes 

(cos0) = {iJo{(pi)\ COS (pi\iJo{(pi)) . (23) 

These equations may be solved for (cos0) and Tc, where the critical temperature can be 
determined by (cos0) — » 0. 

III. RESULTS AND DISCUSSION 

We have considered long-range phase coherence among wires in the bundle in order to see 
whether the phases on the wires are coherent and the bundle as a whole is superconducting or 
not. The self-consistent equation gives rise to a phase diagram exhibiting superconductivity, 
which can be defined as the greatest temperature and field such that cos 6 takes on a non-zero 
value [20]. Here, we assume that the Josephson coupling is independent of a temperature. 



We consider the temperature dependence, Vl — t for p. In order to simplify our calculations, 
we consider the case of the periodic boundary condition. 
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A. No magnetic field 



We consider the following self-consistent equation, substituting B^ff = for the differ- 
ential eq. dsn), 

O 2^ - 4:ZnPl J {cos (f)) COS0i + ^Z^Jpl^l + (cOS0)^) V'n(0i) = En'4)n{.(t)i) ■ (24) 

This equation can be reduced to the standard Mathieu equation [21], using t> = 0/2, y{y) = 

V^n(0./2), 



+ {an - 2qcos2v)yn{v) = 0, 



(25) 



where the characteristic value of the Mathieu equation and q are written as 



4{En-2ZnJpt{l + {cOS(j)y)) 



2,,2mp2 5(1 + (cos 0)2) 



2 2 T/ J.\2"^P0 
iZnPoJ {cos (p) - 



B 



(COS0), 



A 



f9 A 

where we define A = and B = 2znJpl. The eigenvalues are explicitly written as 



En = Aan + B{l + (cos 0)2). 



(26) 



The allowed eigenfunctions are determined by the condition that the wave functions be 
single- valued, i.e., that ipn{(p + 2tx) = ipn{<P)i or equivalently, that ?/„(f + tt) = yn{v). The 
allowed three lowest solutions, up to the order of q'^, are [21 1 



yo{v,q) = i 

'TT 



y2{v,q) 



y-2{v,q) = 



1 — - COS 2f + g 



cos4f 
32 



le) 



2 



cos2t> — g 



sin 2v — q 



cos4t> 



12 
sin4f 



cos6f 19cos2t> 



384 



288 



02 = 4 + 



12 ' 



sin 6f sin 2v 



a_9 = 4 , 

12' 



12 ' V 384 288 

where these are normalized like J^^ ■ipn{(f))d(j) = 1. Thus, the matrix elements for cos6 on 
the corresponding bases, n = 0, 2, and —2, are 



^ ^ 4. ^ 



(cOS0) 



v 



2 V2 

J_ 59 n 

V2 12 ^ 

0-^ 



(27) 



12/ 



From the mapping, we can get the self-consistent condition for the critical temperature 
of the phase ordering in terms of the parameters of the GL equation for sufficient or infinite 




FIG. 1: Temperature dependence of the phase-ordering for mfinite long wires. 



long wires eq. (123|) . which corresponds to the only consideration of the ground state (n = 0) 
in the quantum mechanics problem, and it takes the following form. 



COS( 



A 



2' 



(28) 



The temperature dependence of order parameter obtained by eq. (1251) for infinite long wires 
is shown in Fig. [TJ This figure clearly shows that there is a second order phase transition 
at t = 0.5 because the order parameter continuously becomes zero at the critical point. As 
expected, the critical temperature of the whole wires is lower than the critical temperature 
of a single wire. The transition temperature of phase ordering can be calculated by finding 
the temperature where (cos0) becomes zero. Thus, because with A — > gfj^2%i^f-j ^^'^ ^ ^ 

2z„KA'^(t) 
f5 , 



B 
A 



where a = ^^t^tt-^Sw; the condition becomes 



(cos0) = — > tc = y/c^g'^itc 



Therefore, using g{t) = a/1 — t, this critical temperature becomes 

' 1 + v^' 



(29) 



(30) 



(31) 
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On the other hand, for finite length wires, contributions from excited states in the quan- 
tum mechanics problem need to be considered because the effective temperature is not zero. 
Using up to the order \n\ < 2 for the solution of Mathieu's equation, using eq. (122]) . the 
following self-consistent condition can be obtained, 

-—n= 12l /02^ 

e-/3.//£;o + e-/3.//S2 + e-/3e//i?-2 ' ^^^^ 

where PeffEn = /?e//(^«n + -B(l + (cos0))), but the second term can be canceled. Therefore, 
with the mapping, we can get 

(33) 



i-t 



where we use the following mapping 

" l-t ° eo 1-t A2(0) ^0 l-t ^ ^ 

Using the numerical values according to Tang et al py], Xq ~ 1.4 x 10"^. A plot of T^. versus 
a for several lengths (lOO^o; lOOO^o; 2000^0; and SOOO^o) and infinite length are given in 
Fig. [21 This figure shows that as the length of the wires has increased, the phase critical 
temperature has increased. 
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B. Perpendicular magnetic field 



The critical temperature for the presence of a magnetic field on the wires can be obtained 
by solving the non-Hermitian eq. fl?I]) . 

-h^ e*B,ffh d 



2mpQ ( 



2mc i 



-AZnPoJ (cos 0) COS </)j-+ 2ZnJpoil+ (cOS (f)) ) > ^/'„(0i)= 



Using 'ipn{4>) = eP'"F{v) with and v = 0/2 and 



p = i- 



ch 



again this equation reduces to the standard Mathieu equation: 

d^F(v) 



{2qcos2v)F{v) = -a^F{v) 



(35) 



(36) 



(37) 



where 



a^-p' = A{En-2znJpt{l + {cos(f)y)) 



iZnPoJ {cos (p)- 



B 



(cos0). 



A 



h? A 

The allowed eigenvalues are determined by the boundary condition that ipni'P + 27r) = 
iljn{4>), or equivalently F{v + vr) = exp{—pTT)F{v). Thus we are interested only in the 
Floquet solutions of the Mathieu equation with Floquet exponent u = 2n + ip, where 
n = 0, ±1, ±2, These solutions are explicitly written as 2l| 



F^{v) = coe' 



1-q 



-,2iv 



4(z/ + l) 



(38) 



where Cq is a normalization constant. The eigenvalues are, using g = — 4 cos 6* 



2(z/2 _ r 



(39) 



The allowed three lowest solutions, up to the order of g^, are |21[] 




TT 



COS 2v + p sin 2v 
2(1 +p2) 



, dip - 



2{l+p 



\n\ 4:\3 + ip 1 + ipJJ 2(- 



4 \3 + ip 



+ Aip + 3) ' 



2(3-p2_4ip)' 
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Left wave functions can be obtained from right wave function with ilj^{v,p) = ipniv-, —p)*- 
The self-consistent condition for long wires becomes the following form, 

A 



(COSi 



q 



2(1 +p2)' 

because the matrix elements for cos6 corresponding to n = 0,2, and —2 are, using q 

B 

A 



(40) 



^■(cos^). 



COS( 



v 



1+2? 
1 

1 



3+4jp-p2 2{l+p2) 
1 2 , 

2(l+p2) 3_4ip_p2 y 



(41) 



Again, we can determine the transition temperature of the phase ordering, where (cos 
becomes zero. 

= 0^t, = 



(COS) 



a 



l+p2(tc 



:9 (tc 



The approximation g{t) = vT^^ is again used for this case and then we can get 

t. 

where we define / as 
P 



1 + V^^^' 



(42) 



(43) 



- Jo- 



I'Kh^ A2(0)A,^ol-i 



(44) 



where $o = hc/e* . When / = 0, this solution corresponds to the previous case. A plot of Tc 
versus a for / = 0, 0.2, 0.4, 0.6, and 0.8 is given in Fig. [31 This figure shows that the critical 
temperatures have the minimum values for the interaction between the wires. These values 
can be calculated by 



z„ K > 



The critical /c, which is related to the maximum flux in the wires can be obtained. 



fc 



^a(l-t)2-t2 



(45) 



(46) 



Near the critical temperature of phase ordering, using t = - 5t, this can be written 

fc ^ 72(1 + ,/^)a^/^V6i. (47) 

Fig. m shows that this critical magnetic field fc for a = 0.2, 0.4, 0.6, 0.8, and 1 as a function 
of a temperature. 
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a=1 




0.2 0.4 0.6 

Temperature TIT^ 

FIG. 4: Temperature dependence of the critical field strength, fc for a = 0.2, 0.4, 0.6, 0.8, and 1 
for infinity long wires. 
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IV. SUMMARY 

We have presented a mapping between a one- dimensional GL problem in the presence of 
a vector potential along wires and a two-dimensional quantum mechanics problem with a 
perpendicular magnetic field. Moreover, in the case of weak links between wires, we have 
obtained, using the mean-field approximation, the phase diagrams for the presence of a 
magnetic field and absence of it. 

Next, we discuss the parameters used in this paper. Using the numerical values of the 
various parameters appropriate to those of a single- walled carbon nanotube, which according 
to Tang et al [1] , where superconducting with a relatively high transition temperature T° = 
15K, keT^ = 1.3meV, aoT^ = 6meV, 7 = l.SmeVA, m* = 0.36m^, and ^0 = . " = 
42A, we can obtain the following values for a and f, a = „ ^ f"''^ — ^ and f = 1.7 x 10^^%^. 

' ° ' 8.6xlO~°meV <I>o 

The Josephson coupling energy K is approximated by 2m s'^ where s is the distance between 
nearest wires. If we use rric = m*, K can be written as ° ^2 " ■ Thus, supposing s ~ 5A, K 
is order of 100 ~ 1000 [me V]. Therefore, our values used in the figures are well suited for 
describing real systems. 

We discuss about the use of the GL free energy functional. In principle, this free energy 
functional is applicable only near the critical temperature, T — T° ^ T°. Besides near the 
critical temperature T°, the qualitative description of this functional may not be reasonable, 
although we can employ higher order expansions of the order parameter in the G.L equation. 

We want to comment the effect on the interaction term by a magnetic field. When there 
is a magnetic field, the phase difference needs to be replaced by (pi — (pi^i — / A- dl where 
the integration is between different wires. However, because the direction of vector potential 
is taken in the direction of the wires, there is no contribution from the integral on the 
phase difference. 

In this paper, we only consider the periodic boundary condition for simplification. When 
wires are sufficient long, the effect of the boundary conditions may not change the physical 
properties of the system. However, these boundary conditions may affect the properties of 
the system because of finite length of wires. Moreover, our theory neglects the effects of 
disorder, which plays an important role on balk superconductors. With these degrees of 
freedom, the properties of the system may be changed. Thus, it might be an interest to 
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consider these cases for our future research. 
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